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We investigate the spontaneous emission rate of a two-level quantum emitter near a graphene-coated substrate 
under the influence of an external magnetic field or strain induced pseudo-magnetic field. We demonstrate that 
the application of the magnetic field can substantially increase or decrease the decay rate. We show that a 
suppression as large as 99% in the Purcell factor is achieved even for moderate magnetic fields. The emitter’s 
lifetime is a discontinuous function of |B|, which is a direct consequence of the occurrence of discrete Landau 
levels in graphene. We demonstrate that, in the near-field regime, the magnetic field enables an unprecedented 
control of the decay pathways into which the photon/polariton can be emitted. Our findings strongly suggest 
that a magnetic field could act as an efficient agent for on-demand, active control of light-matter interactions in 
graphene at the quantum level. 


The possibility of tailoring light-matter interactions at a 
quantum level has been a sought-after goal in optics since the 
pioneer work of Purcell!, where it was first shown that the 
environment can strongly modify the spontaneous emission 
(SE) rate of a quantum emitter. To achieve such objective, 
several approaches have been proposed so far. One of them 
is to investigate SE in different system geometries*-!!. Ad¬ 
vances in nanofabrication techniques have not only allowed 
the increase of the spectroscopic resolution of molecules in 
complex environments!*, but have also led to the use of nano¬ 
metric objects, such as antennas and tips, to modify the life¬ 
time, and enhance the fluorescence of single molecules!!!!. 
The presence of metamaterials may also strongly affect quan¬ 
tum emitters’ radiative processes. For instance, the impact of 
negative refraction and of the hyperbolic dispersion on the SE 
have been investigated 17 —. Also, the influence of cloaking 
devices on the SE of atoms has been recently addressed 20 . 

Progress in plasmonics has also allowed for a unprece¬ 
dented control of light-matter interactions at a quantum level. 
When the emitter is located near a plasmonic structure it may 
experience a strong enhancement of the local field. This ef¬ 
fect can be exploited in the development of important appli¬ 
cations in nanoplasmonics!! 23 . However, structures made of 
noble metals are hardly tunable, which unavoidably limit their 
application in photonic devices. To circumvent these lim¬ 
itations, graphene has emerged as an alternative plasmonic 
material due to its extraordinary electronic and optical prop¬ 
erties^—. Indeed, graphene hosts extremely confined plas- 
mons, facilitating strong light-matter interactions!!!!. In ad¬ 
dition, the plasmon spectrum in doped graphene is highly tun¬ 
able through electrical or chemical modification of the charge 
carrier density. Due to these properties, graphene is a promis¬ 
ing material platform for several photonic applications, spe¬ 
cially in the THz frequency range!!. At the quantum level, the 
spatial confinement of surface plasmons in graphene has been 
shown to modify the SE rate!*!!. The electromagnetic (EM) 


field pattern excited by quantum emitters near a graphene 
sheet!! further demonstrates the huge field enhancement due 
to the excitation of surface plasmons. A graphene sheet has 
also been shown to mediate sub- and super-radiance between 
two quantum emitters!!. Recently, the electrical control of 
the relaxation pathways and SE rate in graphene has been ob¬ 
served!!. Despite all these advances, the achieved modifica¬ 
tion in the emitter’s decay rate remains modest so far. Most 
of the proposed schemes consider emitters whose transition 
frequencies are in the optical/near infrared range, usually far 
from graphene’s intraband transitions. As a consequence, the 
effects of graphene on the SE rate are only relevant when the 
emitter is no more than a few dozen nanometers apart. 

Here, we propose an alternative mechanism to actively 
tune the lifetime of a THz quantum emitter near a graphene 
sheet by exploiting its extraordinary magneto-optical proper¬ 
ties. We show that the application of a magnetic field B al¬ 
lows for an unprecedented control of the SE rate for emitter- 
graphene distances in the micrometer range. This is in con¬ 
trast to previous proposals, in which the modification of the 
SE rate was achieved by electrically or chemically altering 
graphene’s doping level. The fact that we consider a low- 
frequency emitter enables us to probe the effects of intraband 
transitions in graphene on the decay rate, which have also 
been unexplored so far. In summary, our key results are (i) 
a striking 99% reduction of the emitter SE rate compared to 
the case where B = 0; (ii) a new distance-scaling law for 
the decay rate that corrects the typical 1/d 4 behavior and is 
valid for a broad range of distances and magnetic fields; (iii) a 
highly non-monotonic behavior of the SE rate as a function of 
|B|, with sharp discontinuities in the regime of low tempera¬ 
tures; and (iv) the possibility of tailoring the decay channels 
into which the photon can be emitted. These findings can be 
physically explained in terms of the interplay among the dif¬ 
ferent EM modes and of electronic intraband transitions be¬ 
tween discrete Landau levels in graphene. 
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I. METHODS 


Let us consider the situation depicted in Fig |T] The half¬ 
space z < 0 is composed of a non-magnetic, isotropic, and 
homogeneous material of permittivity £ s (u>), on top of which 
(z = 0) a flat graphene sheet is placed. The system is under 
influence of a uniform static magnetic field B = Bz. The up¬ 
per medium z > 0 is vacuum and an excited quantum emitter 
is located at a distance d above the interface. 



Figure 1. Quantum emitter at a distance d above a graphene sheet 
on the top of a substrate of permittivity e s (co). The whole system is 
under the influence of a magnetic field B = Bz. 


We consider that the quantum emitter dynamics is 
well described by two of its energy eigenstates (|g) 
and |e)). Within the electric dipole approximation 
and weak-coupling regime, one can show that the SE 
rates I'x and T|| for transition dipole moments perpen¬ 
dicular and parallel to the XY - plane respectively are 
(Appendix A )£r— 
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where To = | d eg | 2 loq/{^'K£qHc 3 ) is the free space SE rate, d eg 
is the emitter’s electric dipole matrix element, ujq = koc is the 
transition frequency, k z = y/k^ — k 2 , and r s,s , r p p are the 
graphene-coated wall polarization preserving reflection co¬ 
efficients. Although the cross-polarization reflection coeffi¬ 
cients r s,p and r p s are non-vanishing in the case of graphene 
under the influence of an uniform static magnetic field, being 
responsible for Faraday and Kerr rotations, they do not con¬ 
tribute to the emitter’s lifetime in the present situation (see 
Appendix A). The diagonal reflection coefficients are given 
by (Appendix B A— 

,, A 2 + A^A T nn A 2 + A L Al 

A 2 + A+A+ ’ A 2 + A+A+ ’ k ’ 

where A 2 = Z^k z k z afj, A± = k z £ s /e o ± k s z + 
k z k s z cr L /(u £ 0 ), = k z ± ki + nou>a T , Z 0 = y/fj^Je o, 

k s z = y/fj, 0 £ s u>$ - k 2 , and k = |fc| = \k x -x + k y y\. 


<jl, <7t and cth are the longitudinal, transverse and Hall 
conductivities of graphene, respectivelly, which are in gen¬ 
eral functions of both frequency and tranverse wavevec- 
tor k. Although the dependence of the material proper¬ 
ties on wavevector may be relevant in the near-field^, we 
have checked that this is not the case for the distances 
we consider. Indeed, the evanescent waves contribution to 
the SE process is suppressed by a e~ 2kd factor, whereas 
non-local effects on graphene’s conductivity become signif¬ 
icant for k > max (^y/eWfh, luq/vf, Here, 

vf — 10 6 m/s and r is a phenomenological relaxation 
time of electrons in graphene. Therefore, provided d 3> 

min (^y/K/eB , vf/ojq, Vft^J we can safely set k = 0 in the 
conductivities, in which case er F = a?. 

We will study the lifetime of quantum emitters in the low 
temperature (fcgT <C /x c ) and frequency (Hcoq -C /r c ) regimes, 
where fj c is the graphene’s chemical potential. As a result, 
graphene’s conductivities can be approximated by their intra¬ 
band terms 
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where Ai ntra = M„ c+ i — M rlc , M n = y/nM\ are the Landau 
energy levels, M 2 = 2heBv F , and n c = int[p, 2 /M 2 ] denotes 
the number of occupied Landau levels. 


II. RESULTS 


Following previous experimental work on SE^, we con¬ 
sider from now on an emitter with a strong transition at u>o = 
4.2 x 10 12 rad/s (~ 0.7 THz). We set r = 0.184 p A \i c = 115 
meV and, inspired by recent experiments on magneto-optical 
effects in graphenA, consider a silicon carbide (SiC) sub¬ 
strate. It is important to clarify that ujq may be a function of 
d (the emitter energy levels can be Lamb-shifted) and of B 
(the levels may also be Zeeman-shifted). However, for the 
purposes of the present work, both effects may be neglected. 
A numerical estimate shows that for the distances considered 
here, the influence of the Lamb shift on the SE rate is unno- 
ticeable, regardless of the value of B. Concerning the Zeeman 
shift, we have checked that although some energy levels may 
be altered in their absolute values, the suppression and en¬ 
hancement factors of the SE rate due to the application of B 
are insensitive to this shift. 

In Fig. [2]we plot the normalized SE rate Fx/Lo as a func¬ 
tion of the distance d between the emitter and the half-space 
for several values of B. For d > 100 //rn the coupling be¬ 
tween the emitter and the graphene-coated wall is mediated 
by propagating modes (k < ko) of the vacuum EM field. 
In this regime of distances the emitter’s lifetime is barely 
affected by B. This behavior results from the fact that in 
the far-field the phase e 2lkzd gives a highly oscillatory inte¬ 
grand in Eq. 0 , except for k z ~ 0. In this case, however. 
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Figure 2. Normalized decay rate Tx/Eo as a function of distance 
d between the emitter and the graphene-SiC half-space for different 
magnetic fields. The inset presents the relative SE rate AEx as a 
function of d for the same values of B. 


r s,s ^ f p,p ^ —1 + O(k~/ko), so that the reflectivity of 
the graphene-coated half-space is almost saturated. Hence B 
hardly affects the reflection coefficients in this regime. A tran¬ 
sition from the oscillating pattern at large distances to a sharp 
growth at small distances takes place for d < 100 //m. In this 
regime of distances the emission is dominated by evanescent 
modes (k > ko) of the vacuum EM field. Interestingly, for 
d < 10 //m changing B strongly affects the lifetime of the 
quantum emitter. A striking suppression of 99% in the Pur¬ 
cell factor, when compared to the case where B = 0 T, occurs 
for 1 jim < d < 10 //m and B > 10 T. Even for smaller 
values of B the Purcell effect is greatly reduced. For exam¬ 
ple, for d = 3 /rm the influence of the graphene-coated wall 
on Tx can be reduced by a factor of 10 for B = 5 T. These 
results are highlighted in the inset of Fig. [2] where we plot 
ATx = [r_|_(<i, B) — Tx(rf, 0)]/rj_(d, 0) as a function of d 
for the same values of B. For even smaller distances an en¬ 
hancement of the SE rate takes place as the magnetic field in¬ 
creases. For clarity this effect is not shown in Fig. [2] although 
it can be noticed in the inset for d < 1 /rm. For instance, for 
B = 5 T and d = 0.2 /rm the SE rate is enhanced by ~ 500%. 

It is also interesting to analyze the distance-scaling law of 
the SE rate for graphene under an external magnetic field. In 
the near-field regime, one can show that (Appendix C) 

1% _ 3eoC 3 Re[cr L ] 1 / |Imox| \ 

r 0 “ u$(e s + £ 0 ) 2 d 4 V w o(e s +£o)<V ’ 

where F(x) is defined in Appendix C, provided Re [ox] -C 
wo(e s + £o)d and Im[e s (wo)] — 0. The validity of this equa¬ 
tion is not restricted to the case when a magnetic field is 
present, rather it is valid whenever correction due to Im ox 
are non-negligible. Equation ({6} explains the results in Fig. 
[2] for a broad range of distances (0.3 /rm < d < 1.4 /rm) 
and magnetic fields (5 T < B < 20 T) with error < 10%. 
The distance scaling-law Ti oc d~ 4 F(do/d) [where do = 
|Imox|/wo(£ s + £ o)] differs from the recently observed result 
Tx oc d ', obtained in the case B = 0--. This difference 
arises due to (low frequency) intraband transitions and losses 
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Figure 3. The decay channel probabilities as a function of (a) d, 
and (b) B for /r c = 115 meV. In (a) solid and dotted curves are for 
B = 5 T and B = 15 T, respectively. In (b) the distance is fixed at 
d = 4 gm. Tx (d, B)/To as a function of B is plotted in (c) d = 200 
nm, and (d) d = 1 /rm. The vertical lines show the position of the 
peak of Tx(d, B)/Tq obtained via Eq. <1 111 . 


in graphene, whose signature is coded in the function F(x ) 
appearing in Eq. (|6j. Indeed, while in the high frequency 
regime (wo r _1 ) graphene’s conductivity is approximately 
a real function, this is not true at the frequency considered 
here (wo ~ r _1 ). However, the Tx oc d~ 4 can be derived 
provided |Im ox| -C wo(e s +£o)d- Since Im[ox] is greatly af¬ 
fected by B, the magnetic field could be exploited to tailor the 
distance ranges where this condition is satisfied, allowing for 
a real time control of the distance-scaling law in the near-field. 
Note that the effects of B on the SE are predominantly related 
to changes in ox. We have verified that 07/ can be neglected 
in Eq. i[3} for the chosen material parameters. In this case, 
the same modifications in the SE rate could be obtained by 
applying a trigonal distortion to graphene, which would gen¬ 
erate a strain induced pseudo-magnetic field, giving rise to the 
formation of Landau levels in graphene’s electronic spectrum, 
while keeping an = 0 due to time-reversal invariancei^s 

To understand such an influence of B on the SE rate in the 
near-field it is necessary to delve a little deeper into the decay 
process itself. The spontaneous decay of a source is often as¬ 
sociated to the emission of radiation to the far-field, but that is 
not necessarily the case. In particular, in the near-field regime 
the emitter decays preferentially into non-radiative channels, 
like surface waves characterized by k > ^/e s / e^k^. For the 
transition frequency we are considering, |Re ox| ~ |Im ox|, 
so that surface magneto-plasmon polaritons^ are strongly 
damped, playing essentially no role in the SE process^. Nev¬ 
ertheless, the so called lossy surface waves (LSW ) 42 ' 55 ' 56 are 
crucial here. These waves correspond to non-radiative pro¬ 
cesses and emerge in the case where the emitter’s energy is 
transferred directly to the substrate or graphene, generally giv¬ 
ing origin to an excitation (e.g. electron-hole pair). Such 
waves are quickly damped, with their energy being usually 
converted into heat 42 —. In the extreme near-field regime, ab¬ 
sorption in the materials governs the SE process and the LSW 
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(k A'd) are usually the main channel into which the emitter 
loses its energy. 

In Fig. 0 we unveil the role played by the different decay 
channels in the emitter’s lifetime. Fig. 0 a) depicts the decay 
probabilities p\, P™’ and p^f^ of energy emission in a prop¬ 
agating (P), totally internal reflected (TIR), or LSW mode, 
respectively, as functions of d for two different values of B. 
These probabilities are given by the ratio between the partial 
decay rates into the aforementioned modes and the total SE 
rate. The partial contribution of propagating, TIR, and LSW 
modes to the SE rate can be respectively well approximated 
by (Appendix D) 


T p : 3 f k ° k 3 Re[r p ' p e 2 V fc o- fc2d l 
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where Tq S p are the reflection coefficients of the graphene-on- 
substrate system in the quasi-static limit [equivalent to take 
c —> oo in Eq. (0|. 

We note in Fig. [3] that changing B can severely affect the 
possible decay channels in the 1 pm < d < 10 pm range, 
essentially swapping the role of the LSW and TIR modes as 
the dominant decay pathway. Indeed, for d = 4 pm we note 
that pY SW drops sharply from 75% to 15% when B changes 
from 5 T to 15 T. On the other hand, p™ (Pj_) increases from 
20% (5%) to 67% (18%). This effect is evinced in Fig. 0b), 
where we plot the decay probabilities as a function of B , for 
d = 4 /an. It is then clear the overall downward (upward) 
trend of p^ sw (p™) as B is increased, with a dominance ex¬ 
change at B ~ 10 T. 

Figures0c) and0d) show 1% (d. B) /l’u as a function of 
B for d = 200 nm and d = 1 pm, respectively, and two 
distinct values of p c . The SE rate presents sharp discontinu¬ 
ities, which are directly linked to the discrete character of the 
Landau levels brought about by the application of B. These 
discontinuities occur whenever a given Landau level energy 
crosses Moreover, there exists a critical magnetic 

field B c = p 2 /{2hev 2 F ) above which the discontinuities are 
no longer present. This is due to the fact that for B > B c all 
positive Landau levels are above p c , so no more crossings can 
occur. Note that the curves merge in the final plateau, regard¬ 
less of the value of p c . For B > B c we have Ai ntra = Mi, 
that does not depend on p c . Hence, provided fcgT -C p c both 
<jl and an are approximately independent of p c for B > B c . 
As a function of B, the SE rate presents a maximum whose 
position depends on both p c and d. This behavior can be un¬ 
derstood recalling that for short distances the SE rate i s 3339 


p o r°° 

“ 2 k 3 J dk P( k ) lm [ r ' P ' P ( k ’ U} o,B)] , (10) 

where p(k) = k 2 e~ 2kd has a maximum at fc™ ax = 1/d. 
In the large k limit Im[r p,p ] presents a peak at fc™ ax ~ 



Figure 4. 3D plot of the relative spontaneous emission AFj (p c , B) 
as a function of both p c and B for d = 2 pm. 


£oWo[e s /eo + 1 ]/|ctl|. Since \<jl{uq,B)\ decreases with B 
(for B > 1 T in our case) we note that fc™ ax moves to high 
values of k as B increases. Therefore, for a fixed emitter- 
graphene separation, the overlap between p{k) and Im[r p ’ p ] 
grows with B until /c™ x ~ fc” ax . After that, this overlap 
diminishes and so does r_L, which explains the behavior of 
the SE rate in Fig. 0 The value of the magnetic field B m 
that maximizes r_i_ can be estimated by setting /c™ ax = fc™ ax . 
This leads to 

I°'l( w 0 j Pc, B m )\ — £o w o< 7 [e s ( w o)/ £ o + 1 ] ■ (H) 

The accuracy of this equation is clearly seen in Figs. 0c)-(d) 
where we show B m calculated through Eq. CD for p c = 115 
meV and p c = 150 meV. 

Similar results hold for T || in the near-field regime. In¬ 
deed, for d <C 27r/fco the contribution of r s,s to T|| is neg¬ 
ligible and the approximation k z ~ ik is valid. Hence, 
apart from a factor 1/2, F|| can also be written as in Eq. 
( flOl) (see Appendix C)^. In Fig. 0we plot AT^(p c ,B) = 
[T\\(p c ,B) - r||(// c , 0)]/r||(/x c , 0) as a function of both p c 
and B for d = 2 pm. In this case, the reduction in the Purcell 
factor in the pm range can be as high as 98%, when com¬ 
pared to the case B = 0. Figure 0 corroborates our conclu¬ 
sions that an astounding control on the radiative properties of 
quantum emitters can be achieved via magneto-optical prop¬ 
erties in graphene. Moreover, Fig. 0 reveals that the SE rate 
can be modified by keeping B constant while changing p c , 
which could be implemented by applying a gate voltage on 
graphen a 2(>3 ° . 


III. CONCLUSION 

In conclusion, we have shown that the application of a mag¬ 
netic field allows for a great control over the Purcell effect 
and decay pathways of quantum emitters near graphene. Alto¬ 
gether, our findings demonstrate the viability of actively dic¬ 
tating optical energy transfer processes with magnetic fields 
or strain. By demonstrating that these results are within the 
reach of state-of-the-art experiments on quantum emission in 
the THz range, we expect that they may find further applica¬ 
tions in quantum photonics, and may even serve to probe other 
light-matter phenomena. 
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APPENDIX A: SPONTANEOUS EMISSION NEAR AN 
ANISOTROPIC INTERFACE 


The spontaneous emission decay rate of a two level emitter 
can be written in terms of the EM dyadic Green’s function as 

2w 2 

r = ^^2 Im [ d ge • G(ro, ro; Wo) ■ d ge ] , (12) 

where r 0 = (0,0, d) is the position of the quantum emitter and 
G(r, r'; w) is the EM dyadic Green’s function, which allows 
one to write the electric field at position r that is generated by 
a point dipole at position r', oscillating with frequency w as 

E(r;w) = /xw 2 G(r,r';w) • d(w), (13) 

where // is the permeability of the medium where the dipole 
is embedded. In vacuum the dyadic Green’s function satisfies 
the inhomogeneous Helmholtz equation 

uj 2 

V x V x G(r, r'; w) - -^Gfr, r'; w) = M(r - r'), (14) 
& 

with I being the unit dyad. 

We are interested in the case where the emitter is located at 
a distance d above a semi-infinite homogeneous medium with 
flat surface at z = 0, where the graphene layer lies. Since the 
dyadic Green’s function obeys the same boundary conditions 
as the electric field, we can write it for z, z' > 0 as 

G(r,r';w 0 ) = G (0) (r, r 7 ; w 0 ) + G (r) (r, r'; w 0 ), (15) 

where G (0) (r, r'; wo) is the free space Green’s function and 
G w (r,r> 0 ) is the reflected one. For z' > 0 and z < 0, 
the dyadic Green’s function can be written as a transmitted 
Green’s function, G (t) (r, r'; uiq). Each of these Green’s func¬ 
tions can be conveniently expressed in terms of its spatial 2D- 
Fourier transform G (0/r/t) (fc, z,z';u>) as^i 

G (0/r/t) (r,r';w)= / -^re ik < x - x,) & 0,llt) {k, z, z'; w), (16) 

J (2 tt ) 2 


where k = k x x. + k y y and x = xx + yy. The free space 
Green’s function is given by 

G { 0 ) (k,z,z'-uj) = -^-e lk ’‘ lz ~ z ' 1 (ep ® 6p + ef O e±) ,(17) 
with k z defined as 


k z 


yJk,Q — k 2 , k < k 0 
iyjk 2 — fcg , k > ko 1 


(18) 


and we have introduced the polarization vectors for s- and p- 
polarized waves (the + and — signs correspond to z > z' and 
z < z ', respectively) 

± ky±-k x y ± fc , k z k x k + ky y 

e s = - J —r -, e p=7T z= Fl- T. ” • ( 19 > 

k ^ ko ko k 


Note that these vectors are orthogonal, but they are normalized 
only for propagating modes (k < ko). 

The reflected Green’s function can be written as 


G (r) (fc, z, z'\w) 


1 ik z (z+z') 

2 k z C 


Y. ri,i 4 ® Cj r . (20) 

i,j=s,p 


where r 1 ' 1 are the reflection coefficients for an incoming j- 
polarized wave that is reflected as an i-polarized wave. Simi¬ 
larly, the transmitted Green’s function is given by 


G (t) (fc, 


= 2T e 


—ik s z z gik z z 


' ,t 


e i 


( 21 ) 


where t l,i are the transmission coefficients (incoming j- 
polarized wave, transmitted i-polarized wave) and ef t are the 
polarization vectors in the substrate, given by Eq. 03 af¬ 
ter replacing ko by koyJe s /eo and k z by k z . The reflection 
and transmission coefficients are obtained by imposing the 
usual boundary conditions on the EM field at 2 = 0 and by 
modelling graphene as a two-dimensional current distribution 
(see Appendix B). 

The evaluation of the SE rate requires the evaluation of the 
dyadic Green’s function at the coincidence r' = r = ro. In 
this case the integration over the momentum angular variable 
in Eq. can be easily performed. The only nonzero com¬ 
ponents of G (0) (ro, ro; wo) are the diagonal ones. The contri¬ 
bution of G (r) (ro, ro; wo) to the SE rate presents polarization 
preserving terms (which involve e+ <g> e~ and e+ ® e“) and 
cross-polarization terms (which involve e+Ge“ and e+ 0f p ). 
After performing the angular integration, the polarization pre¬ 
serving terms only select the diagonal terms of d* c 0 d ge . 
The cross polarization terms select the d* c x d gejy — d* e y d g e,x 
components of the dipole matrix elements. As the transition 
dipole matrix elements of a two-level system can be made 
real by a proper choice of the relative phase between |g) 
and |e), the cross polarization terms do not contribute to the 
SE rate. Therefore, only the reflection coefficients r P:P and 
r s,s give a non-vanishing contribution to the SE process even 
though cross polarization coefficients r s,p and ?’ p,s are non 
zero. Similar conclusions hold even in the case of a semi- 
infitite anisotropic substrate. By plugging Eqs. into 

m it is straightforward to obtain Eqs. ()T| and (0). 
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APPENDIX B: FRESNEL’S COEFFICIENTS FOR AN 
INTERFACE COATED WITH A 2D CONDUCTIVE FILM IN 
THE PRESENCE OF AN APPLIED MAGNETIC FIELD 


Let us consider that an incoming arbitrarily polarized EM 
wave propagating in a dielectric medium with permittivity e\ 
and permeability fii, impinges on the flat interface with a sec¬ 
ond homogeneous medium, with permittivity e 2 and perme¬ 
ability fj, 2, occupying the half-space z < 0 coated by a 2D 
conductive film. For an impinging electromagnetic wave with 
frequency uj and in-plane wavevector k, the electric and mag¬ 
netic fields can be expressed as 


E/ 

H 7 


ERi 

1 r 


Zi L 




0 —ik Zt iz i(k-x-u)t) 


ER - £}€+! 


g—ik Z) iZgi(k-x—ujt) 


( 22 ) 

(23) 


where Ej, E p are the transverse electric and transverse mag¬ 
netic incoming amplitudes, respectively. n and n are 
given by Eqs. ® and ([9]) replacing k 0 by k n = Wy/e n \x n 
with n = 1, 2. Z n = y 7 /i n /e„ is the impedance of medium 
n. Similarly, the reflected and transmitted fields are written as 


and 


E_r = 

H = 

Et = 


E R e s,l + E R e p,l 

1 [ El e- - E f 


gik z ,iz gi(k-x—ujt) 


s,l p,l 


gik Zt iz gi(k-x-ut) 


E s T e+ 2 + E v t e+ 2 


g— ik Zj 2 Z gi(k-x—cot) 


H r =^ \e v t e+ 2 - E s t e+ 

Zo l ’ 


e ~ik z ,2 z e i(k-x—ujt) 


(24) 

(25) 

(26) 
(27) 


We should determine the reflected E R p} and transmitted 
E s T ip) amplitudes in order to calculate the reflection and trans¬ 
mission coefficients 


r 1,J = 


Ek 

e\ 


and f’ J = 


E'rp 

M 


(i,j) = (s,p). 


(28) 


where ey = [k x ± + k v y)/ |fe| and ex = ( k v ± - k x y)/ |fe|. 
o’l (&t) is the longitudinal (transverse) conductivity, <jh is 
the Hall conductivity and is only nonzero in anisotropic 
materials such as black phosphorus^- In case of graphene 
we have = 0, but in order to keep the discussion as 

general as possible and due to the rising interest in black 
phosphorus we will allow for a finite a Using Eqs. 
d22l) -(l27l) into Eq. ( l29l ) and Eq. ( l30l ) one can demon¬ 
strate that the reflected and transmitted amplitudes satisfy 
the following equations 


e s i + e s r = e 


- / Ti 




- Z 2 . 

+ ( a 7y m - <r H ) E t- 


(32) 

(33) 


(34) 


(35) 


Considering separately the cases of s and p incident polar¬ 
ization one can decouple previous equations and show that 
Fresnel’s coefficients in the presence of an external magnetic 
field are given as 


p p + A 2 ss A T ARA 2 

r ’ _ AR + A 2 ’ r ’ - ax + a 2 ’ 
■p.p _ Z 2 e 2 2fe Z| iA^ _ H2 2k Zy iA+ 

Zi£ 0 AR+A 2 ’ /i 0 AR + A 2 ’ 

_b,P ,s,p + 

z 1 /x 2 AR + A 2 ’ 
fcifc z ,2 zl ^^ 2 k z ^k za {(j s y y m -cj H ) 

k 2 k Zi i Zi /Xq AgAg+A 2 


(36) 

(37) 

(38) 

(39) 


with 


The reflected and transmitted amplitudes are obtained by 
solving Maxwell’s equations and imposing the appropriate 
boundary conditions on the interface at z = 0. Taking into 
account the presence of a 2D conductive film at the z = 0, the 
boundary conditions that must be satisfied by the EM field are 

z x [Et — Efl — E/] = 0 , (29) 

zx[H T -H fi - H/] = J 2D = <r • E t , (30) 

where J 2 d is a 2D current density that is induced on the con¬ 
ductive field, and <x is the 2D conductivity tensor of the film^I. 
In the most general case (a 2D homogeneous anisotropic ma¬ 
terial in the presence of a magnetic field) the conductivity ten¬ 
sor can be written as 


A± = (k z ,\£-2 ± k z , 2 £i + k z ^k Zt2 (j L /uj) /e 0 , (40) 

A± = {k Zt 2 Hi ± k z> i/x 2 + w/xi H 2 <Jt) /mo, (41) 

A 2 = Zlmn 2 k z ^k z ,2 [oh - (v7 v m ) 2 ] /$■ (42) 

For graphene a,).)) 171 = 0 and, in the case where medium 1 is 
vacuum (e\ = eg, /xi = ^lq) and medium 2 is non-magnetic 
( fi 2 = /io), the reflection coefficients reduce to the ones given 
in the main text, Eq. 0. 

APPENDIX C: DISTANCE-SCALING LAW IN THE NEAR 
FIELD FOR TERAHERTZ EMITTERS 


er = oxen <g) e|| + ar&± ® ex 
+ cr H (e± ® e|| - ey <g» ex) 

+ ^(ei ® ey + ey ® ex), 


(31) In the near field, the main contribution for the SE rate 
in Eqs. 0 and 0 comes from large in-plane wavevectors 
k k{). In this case the quasi-static approximation holds 
(c -A 00) and k z and k s z can be well approximated by ik. 























7 


Besides, the Hall conductivity gives a negligible contribution 
to the quasi-static reflection coefficients so that we can set 
an — 0. Within these approximations the dominant terms 
in T and I’p originate from the polarization preserving trans¬ 
verse magnetic reflection coefficient and can be cast as 


F i 3 r + °° k 2 

■^2 1 ^ Imra ’ (43) 

ih <«> 

where 

p,p _ i (g« ~ £ o) ~ k<JL ,,n 

QS i (e 0 + £s) w 0 - ka L ' 


In a regime where |Re[cri]| -C (e 0 + e s )ojo d the imaginary 
part of rgg can be approximated by 


Im[; 


,P'Pi 

QS J 


2£ 0 uj 0 kRe[a L } 


{(g> + e 0 )wo - fcIm[cr L ]} 


(46) 


Substituting Eq. (l46t into Eqs. (l43l) and ( 1441 ) one can show 
that the Purcell factor in the near-field regime is given by 


r_L _ 3g 0 c 3 Re[g L ] I f |Im[g L ]| \ 

To _ (e 0 4- £ s ) 2 Wq d 4 \(e 0 +e s )oJod) ’ 

r || 3g 0 c 3 Re[g L ] 1 / lIm[g L ]l \ 

To 2 (£q + £ s ) 2 Wq d 4 \(£o + E s )u}q(1 J 

with the function F(x) defined as 


also to the power emitted by a classical oscillating dipole, 
d(t) = de~ lUot + d*e“°*. The classical power emitted by 
such dipole is related to the SE rate of a two-level quantum 
emitter through P = huJo T 37 i 38 , provided we choose for d the 
transition dipole moment of the quantum emitter. The prob¬ 
ability of decaying into a radiative or non-radiative channel 
can be obtained by computing the fraction of the power that is 
emitted by the classical dipole to the far field and the one that 
is dissipated into the materials, respectively. 

The average power emitted by the classical dipole that 
reaches the far field (radiative processes) can be expressed as 

/*27T /*7T 

P Ta .d = lim / d(j) dds\n(6) r 2 r ■ (S(r)), (50) 

J 0 J 0 

where the Poynting vector in the far field is given by 

S(r,t) = ^E(r,t) • E(r,f)f, (51) 

and (...) denotes time average over one oscillation period. 
Here, Z = \J ///.£ is the impedance of the medium. Us¬ 
ing Eqs. <IT3l) and (fSTb . the time averaged Poynting vector 
can be cast as 

(S(r, t)) = 2 ^° |G(r, r 0 ; w 0 ) • d ge | 2 . (52) 

In the limit fc 0 | r — r 0 | 1, the Green function G(r, r 0 ; wq) 

can be evaluated from Eq. (fl6t using the stationary phase 
method (where the fast oscillating phase is given by k ■ x + 
ikz^\z\). The obtained result is 


F(x) 



y 3 e~ 2v 
(1 + yx) 2 


(49) 


APPENDIX D: DECAY CHANNEL PROBABILITIES 


In order to determine the different decay pathways prob¬ 
abilities one must study how the total power emitted is dis¬ 
tributed into the different channels. Two processes can be 
distinguished: (i) radiative decay, which involves the emis¬ 
sion of a photon that can be detected by a far away detec¬ 
tor; [ii) non-radiative decay, where the emitted power does 
not reach the far-field, but is instead absorbed by graphene or 
substrate and creates a material excitation. In order to com¬ 
pute each channel contribution to the decay process we use 
the fact that the total SE rate given by Eq. (IT2l> corresponds 


ml \ ik n | Z | ik n (i—m f X i J A 

G(r, r 0 ;coo) — e V r JG^-k n ,z,d;u} 0 J , (53) 

with k n = ko for z > 0 and k n = k Q ^/e s /e o for z < 0 . 

Using this result together with equations (l20l >. (fTTI i. (l50l ). and 
(l52l) . one can put the total power emitted into the far field by 
the dipole as -1 


p™d = pz+pi°r, ( 54 ) 

where P'Z and T J r < ( i " 1 wn are the average powers emitted into the 
regions z > 0 and z < 0, respectively. Splitting each term into 
contributions from dipole components that are perpendicular 
and parallel to the XY -plane, performing the integral over (1> 
and changing the variable of integration 9 according to k = 
kg sin 9 for and k = ko \Je s /eo sin 9 for we can 

express the average powers as 
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, 1 3 f k ° k 3 


rad,_L 

Po 

pup 


2 4 


dk- 


/o IM 


=i +i r°dk k 


k 0 

'e _ 

o ko\k z 


|2Re [ 


2Re 


e 2ife^d r p,p1 


|r p,p | 2 + |r s,p | 2 


}■ 


| k z 


k 2 

«,g 


| k z 


| k z 


k 2 


k 2 


Prtd™ 3 k 3 


3 r 

4 Jo 


dk 


„— 2 dlmk z ^ 


kg | k z | 


Po 


pdown 3 fko y fcjT 0 ^ k 


El{l' P ’ P l 2 + l f ’ t l 2 }' 


dk 


- 2 dlmk z I 


\K 


k 0 | k z 


'iri { |f ,B l 2 + %r l tP,P l 2 + l* p ’f + %r l* s ’ p l 

/c 0 A : 0 


i&z 


(55) 
, (56) 

(57) 

(58) 


where Pq = WQ|d ge | 2 /( 37 reoc 3 ) = HujqTq is the total power 
emitted in free space. For the power emitted into the z < 0 
region, there are two different contributions: (i) 7t/2 < 6 < 
arcsin(^/e s /e 0 ), or ko < k < fcoyeT/ip, which is usually 
referred to as forbidden light region 31 - and corresponds to an 
inverted total internal reflection process, in which a decaying 
wave that is emitted by the dipole is transmitted as a prop¬ 
agating wave once it reaches the interface (we refer to this 
contribution as P r ad™’ f )i (**) arcsin(-\/£s/£o) < 9 < n, or 
k < ko, which corresponds to the emission of propagating 
waves (we refer to this contribution as P*^ wn,a ). By subtract¬ 
ing the power emitted via the radiative processes (P ra d) from 
the total dissipated power (Ptotai = fouoT), we obtain the 

I 


power dissipated via non-radiative processes (P non -rad)- The 
contributions of the perpendicular em parallel components of 
the electric dipole to non-radiative power can be written in 
terms of absorption coefficients as 


Pnon— rad,L 
Pnon— rad, || 

Po 


f- d~oo r.3„ —2dlmfc, 



where the absorption coefficients are given by 



1- |r p ’ p 

2 + |r s ’ p | 2 + ^| (|f p ’ p | 2 + |f s ’ p | 2 ) 

, k < k 0 


*4p — < 

21m [r p,p ] 

-jg^P + K-'T) 

, k 0 < k < ko^/ss/eo 

(61) 


21m [r p,p ] 


, koyjes/e 0 < k 



and _A S is obtained from Eq. (16Tb by swapping sop. 

Note that for k > ko \fe s /e o the expressions for the non- 
radiative emitted power, Eqs. < l59t and d60t . coincide with the 
expressions for the total SE rate, Eqs. ([]]) and (0. Hence, we 
can interpret the integration region k > ko \Je s /Eq in Eqs. 0 
and 0 as being a contribution to the SE rate exclusively 
due to non-radiative processes. The region k < ko \J e h j Eq 
also contributes to the non-radiative decay, as can be seen in 
Fig. [5] where the integrands of Eqs. ( 1591 ) and (f60l> as a func¬ 
tion of k are plotted. This is only expected as propagating 
waves emitted by the dipole can also be absorbed and dissi¬ 
pated by the graphene layer. It should be mentioned, how¬ 
ever, that the contribution of wavevectors k < ko\fs7feo to 
the non-radiative SE decay is negligible when compared to 
the contribution coming from k > koyJe 8 /eo (see Fig. 0. 
Therefore, the non-radiative decay due to LSW can be well 
approximated by Eq. 0. In the same way, we can approxi¬ 


mate the contribution to the SE rate from k < ko \Ze s /£0 in 
Eqs. 0 and 0 as being exclusively owing to radiative pro¬ 
cesses. As such we can approximate Prad™ (TIR modes) by 

Eq. © and P r d d™’ a + P ra P _l (propagating modes) by Eq. 0. 
Approximations ©-© were tested numerically against the 
exact results and the differences were found to be negligible. 

Finally, we notice that the power emitted by the quantum 
emitter that is absorbed by graphene due to Joule heating can 
be written as 

P g = 2/i 2 Wo Re / ^ 

S J (2t r) 2 

xdgg-t&i (fc, 0,0; w 0 )-er (fc,w 0 )-G (fc, 0,0; w 0 )-d se . (62) 

In the case when Im[£ s ] = 0, non-radiative decay is exclu¬ 
sively due to graphene and it is possible to show that the power 
absorbed by graphene can be written as Eqs. ( l59l > and (l60l) . 
with absorption coefficients given by Eq. ( I 6 TI 1 . 
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Figure 5. Integrands of Eqs. ( 159b and ( 160b for the quantum emitter's 
dissipated power by non-radiative processes. On the left we plot the 
the integrand of Eq. ( 159b . on the right we plot the integrand of ( 160b . 
splitting it into the individual contributions from the absorption coef¬ 
ficients A p and As. The vertical dashed lines mark the points k = ko 
and k = ko^/e a / eq. The values of B = 5 T and p c = 115 meV 
were used. 
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